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The decoration of graphene samples with adatoms or nanoparticles leads to the enhancement 
of spin-orbit interactions as well as to the introduction of symmetry breaking effects that could 
have drastic effects on spin and electronic transport phenomena. We present an analysis based on 
symmetry considerations and examine the impact on the scattering matrix for graphene systems 
containing defects that enhance spin-orbit interactions, while conserving the electronic total angu¬ 
lar momentum. We show that the appearance and dominance of skew scattering, and the related 
observation of valley and/or spin Hall effect in decorated graphene samples, suggests the set of 
symmetries that adatom perturbations should satisfy. We further show that detailed measurements 
of the transport and elastic times as function of carrier concentration make it possible to not only 
extract the strength of the spin-orbit interaction, as suggested before, but also obtain the amplitude 
of the symmetry breaking terms introduced. To examine how different terms would affect mea¬ 
surements, we also present calculations for typical random distributions of impurities with different 
perturbations, illustrating the detailed energy dependence of different observables. 

PACS numbers: 72.10.Fk, 75.76.+j, 72.80.Vp, 03.65.Pm 


I. INTRODUCTION. 

The low energy Dirac fermions residing near the cor¬ 
ners of the Brillouin zone on the honeycomb lattice of 
graphene 1 have inspired the theoretical prediction and 
subsequent experimental realization of a variety of new 
two-dimensional materials 2 3 and novel electronic phases 
of matter in two and three dimensions. 4 ® The ability to 
manipulate and control the electronic and spin proper¬ 
ties of graphene and related materials by decoration)^® 
stacking) 12 * 13 ^ intercalation 14 18 and distortion 19 27 may 
also lead to important technological advances. 

The helical nature of the massless Dirac particles in 
pristine graphene gives rise to interesting phenomena, 
such as Klein tunneling,® and the observation of weak 
anti-localization in this material!- 29 However, it has been 
shown that the transport properties of graphene strongly 
depend on the nature, symmetries and concentration of 
the impurities present in typical samples pUESl Deforma¬ 
tions from impurities may belong to different symmetry 
classes. As such, they could either enhance or weaken the 
signatures expected from the chiral nature of the scatter¬ 
ed of graphene and drastically modify the electronic and 
spin transport properties. 

An interesting example of interactions protected by 
lattice symmetries is provided by the appearance of a 
spin orbit interaction (SOI) that preserves all symmetries 
of the honeycomb lattice^ and leads to the formation of 
topologically protected edge states-anticipating the exis¬ 
tence of topological insulating materials Pin contrast, an 
electric field perpendicular to a graphene sample leads to 
inversion (z -A —z) symmetry brea king , allowing the ap¬ 
pearance of a Bychkov-Rashba SOlP® The small atomic 
number of carbon makes an intrinsic SOI small (a few 
/ieV),® while the Rashba SOI has stren gth of hundreds 
of /ieV for typical values of electric fields .® 

Methods that rely on the artificial functionalization 


of graphene could lead to sizable enhancement of SOIs. 
Successful examples of these methods include the inter¬ 
calation of gold atoms between grap hene and the Ni sub¬ 
strate p3*®l wea k hydrogenation, 8 9 and the deposition of 
large copper or gold nanoparticlesF 11 Other approaches 
are also promis ing, inc luding externally induced curva¬ 
ture of samples[12122 \2M decoration with heavy adatoms 
such as thallium and indium, 43 44 and the deposition of 
graphene on the surface of a strong topological mate¬ 
rial. 17 These methods have achieved or project the en¬ 
hancement of SOI of up to several tens of meV, lead¬ 
ing to clear measurable effects on the electron dynam¬ 
ic jmnuniisHiH] 

The methods used to enhance SOIs in graphene require 
the presence of adatoms or local defects/deformations 
that accommodate on the graphene lattice in different 
ways. These perturbations introduce additional effects 
from the lattice symmetries they may break. The en¬ 
hancement of intrinsic and/or Rashba SOIs on the dy¬ 
namics of electrons in graphene has been analyzed for 
the most part in the framework of a single-valley Dirac 
description! 8 * 10 * 45 ^® In that limit-valid for potential scat¬ 
tered with large characteristic spatial variation com¬ 
pared to the lattice constant-the scattering preserves he- 
licity, it is predominantly forward, and gives rise to Klein- 
tunneling, as mentioned aboveP- 5 36 45 * If the scattering 
centers lead to enhancement of Rashba and/or intrinsic 
SOIs, the electronic scattering has been shown to become 
increasingly isotropic.® This scattering could also lead to 
spin Hall effect via enhanced skew scattering close to the 
resonant regime P® 

Here, going beyond the single-valley Dirac 
treatment,® we study the scattering properties of 
Dirac fermions from adatoms or imperfections on 
graphene, as they may lead to the reduction of lattice 
symmetries. Through a symmetry analysis and consid¬ 
ering intra- and intervalley processes, we are able to 
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identify the constants of motion for different situations. 
We further examine the constraints imposed by these 
conserved quantities on the scattering matrix and the 
observable consequences in transport experiments. 


II. SYMMETRY OF PERTURBATIONS IN 
GRAPHENE 


A. Spinless scattering 


We focus in particular on the ratio of transport to 
elastic times, £ = r tr /r e , and its dependence on carrier 
concentration, which was previously studied for systems 
where SOIs were locally enhanced, and no additional lo¬ 
cal symmetry breaking effects were considered.^ We find 
that this ratio depends strongly on the microscopic po¬ 
tential perturbations, which make £ deviate from its char¬ 
acteristic low-energy value of 2 P 5 This deviation depends 
on the carrier concentration and strength of the sym¬ 
metry breaking, as well as on the local potential shift 
produced by the impurities. As backscattering becomes 
stronger, the ratio takes a value closer to 1 (the limit for 
massive Schrodinger particles). We further study the ef¬ 
fects symmetry breaking perturbations may produce on 
the detection of Rashba and intrinsic SOIs. The de¬ 
pendence of £ on carrier concentration in the presence 
of Rashba SOI is shown to be qualitatively different for 
different symmetry breaking terms in the perturbation. 
Systematic experimental studies of £ in given graphene 
systems would allow the identification and simultaneous 
evaluation of the relative strengths of the different per¬ 
turbations. Finally, from the symmetries of the scatter¬ 
ing matrix, we show which systems may display non-zero 
valley and/or spin transport skewness (skew scattering). 
In particular, we find that disorder that is not invariant 
under “effective time reversal” leads to a reduction of 
spin skew scattering, due to the enhancement of valley 
skew scattering; this allows one to identify what kind of 
decorated system will display skew scattering and may 
result in the appearance of valley and/or spin Hall effect. 


In order to analyze the effects of symmetry breaking 
on the dynamics of Dirac fermions in graphene, we first 
consider the Hamiltonian 

H = Ho+ V T , (1) 

where 

H 0 = vpa • k , (2) 

(h = 1 and k = (k x ,k y )) is the pristine graphene Hamil¬ 
tonian, represented in the four-dimensional chiral basis 
V’ = {'^a,k^b,k^b,k'^a,k') T , with 

a^r 3 0a j «( % -%) ’ (3) 

(i = 1,2,3), where cq and r* are a set of Pauli matrices 
that act on the sublattice index (A, B ) and the valley 
index (A, A'), respectively. It is also useful to define 

/ = <TO0TO= ( < O <7o ) ’ (4) 

7 5 =T3®<7 0 = (^° , (5) 

and 

/3 = r 1 ®<T 0 =(/ o°) • (6) 


We further illustrate the impact of the symmetry con¬ 
siderations described above for typical graphene samples. 
Throughout the discussion of the different perturbations, 
we present numerical results for the anticipated depen¬ 
dence of £ and skewness to different effects, by analyzing 
the results for a random impurity distribution. This al¬ 
lows us to gain insights into the contribution of different 
effects to various measurable quant it it es. 


The remainder of the paper considers the symmetries 
of different perturbation terms, especially their invari¬ 
ance under time reversal, and the consequences of these 
symmetries on the corresponding scattering matrix in the 
problem. In Sec. m we consider the spinless case, and 
study the scattering matrix and scattering amplitudes, 
as well as their impact in transport and elastic times. In 
Sec. IIB we include the spin degree of freedom and asso¬ 
ciated spin-orbit interaction terms, while Sec. |II C| studies 
the consequences on skew scattering. A general discus¬ 
sion of the results in Sec. Ill is followed by conclusions 
in Sec. IIV1 


As in the case of the three dimensional Dirac equa¬ 
tion EMU these matrices satisfy 


, Oij ) — i 

(7) 

[<*», 7 5 ] = 0 , 

(8) 

IT 

ii 

o 

(9) 

{ 7 5 ,/ 3 } = 0 . 

(10) 


1. Symmetry considerations. 

Here we defin^ Vt, as the set of time reversal invariant 
perturbations ! MMmU 

Wt, 7 ~] = 0 , ( 11 ) 


where 


T = / 3 7 5 aiC 


( 12 ) 
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is the time reversal operator, T 2 = L and C is the anti- 
linear complex conjugation operator. 45 52 53 Additionally, 
the helicity operator 

h = y 5 a • k/k , (13) 


(with k = |fe|) indicates the pseudo-spin projection in the 


direction of momentum,with 

H 0 ,h =0, 

(14) 

and 

{T,h} = 0. 

(15) 


The conservation of total angular momentum in a sys¬ 
tem in which perturbations Vr are present requires 

[H,J Z \= 0, (16) 


where 

J z = -id e + . (17) 

In systems where the spatial dependence of the pertur¬ 
bations has circular symmetry, such as the simple form 
H = Ho + VrO(R — r), there exist a subset of time rever¬ 
sal invariant interactions that conserve the total angular 
momentum, 


Vt,j C Vt • 


(18) 


This set can be_ shown to be composed of the ele¬ 
ments 


V T ,j = {vl, sa 3 , i/3e n5 q , (19) 

where v, s, t and <p are real parameters. Notice that 
with the exception of r>/, the Vt,j terms anticommute 
with aq and <a 2 , so that they generate gaps in the spec¬ 
trum and are referred to as “massive terms” The v 
term describes a constant local potential shift caused by 
the adatoms, while the s term describes a staggered po¬ 
tential effect some perturbations may add, where the A 
and B sites are affected differently. The t terms de¬ 
scribe a hopping modulation amplitude caused by the 
adatoms/defects, where describes the angle of dimer¬ 
ization of the resulting deformation, as in a Kekule pat- 
temPmiH 

One can also define the subset 

Vr,j,h C Vt,j , (20) 

for perturbations that in addition to commuting with T 
and J 2 , also commute with the helicity operator h, and 
find that 


V T ,j,h = {vI,tPe in > t +}. (21) 

In these cases the eigenstates of the Hamiltonian, 


can be identified by corresponding quantum numbers, 
such that 


H±,j = ±ip±,j 


(as h 2 = 1), and 


(23) 


*o> 

-H 

•Os 

II 

-H 

(24) 

where j is half-integer. One can also find an additional 
operator 

Ti = , 

(25) 

such that Tj 2 = — 1, 


[T 1 ,H] = [T 1 ,h] = 0, 

(26) 

while 


{Ti,J z }=0. 

(27) 

7i can be seen as a combination of a chiral rotation 
and the effective time reversal operator in each Dirac 
point JS3tm This operator together with the time reversal 
operator T, relate the states with angular momentum j 
and — j, which can be chosen such that 

Ti ip±,j = ±(-i y^±,-j, 

(28) 

and 


Tip±,j = t(-i y + h T ,-j 

(29) 


(see appendix [A|. 


2. Symmetry constraints on scattering and electronic 
transport. 

Conservation of total angular momentum allows one 
to use the j block-diagonal character of the scattering 
matrix. The lack of valley mixing in the area outside 
the potential perturbation region, r > R, enables the 
determination of the scattering amplitudes between the 
+ and — (or correspondingly K and K') states for each 
j block. The components of each 4x4 block (Sj) are 
given by 

Sj,TT' Sj,TT'&0 1 (^ 0 ) 

where r, r' = d=. Unitarity of S requires 

SjS] = I , (31) 

and symmetry under operators 7i and T determine rela¬ 
tions between the j and — j blocks, such that 

I'TiSjTfV = Slj (32) 

and 


H = Ho + V T ,j,h 


( 22 ) 


7 S T5 i r- 1 7 5 = 


(33) 
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conse- 


(see appendix [A]) . These relations have important 
quences on the possible scattering processes (see Eq. A13 
in appendix), which allow us to obtain different scatter¬ 
ing functions. 

For example, the differential cross section a(6) contains 
information on the angular distribution of the scattering. 
From the relations in Eq. |A13[ one can determine the 
angular dependence of the differential cross section in 
terms of the j-scattering amplitudes. For a circularly 
symmetric region where 


H = H 0 + V T , 


JM 


for r < R, and 


H = H. 


for r > R, we have 

<J TT (0) = 

cr Tf (6) = 


2 

7 xk 

2 

7 rk 


f tt ij cos j 9 

OO 

F frr,j sin jO 

i>2 


(34) 


(35) 


(36) 


<T ( 0 ) =^2 a TT'(0) , 

r,r' 

where f TT ',j are the scattering amplitudes defined in ap¬ 
pendix [Bj and f = —r. Notice that the differential cross 
sections for valley-preserving processes, <j tt (0 ), display a 
maximum at 6 = 0, and vanish at 0 = n (which signals 
the absence of backscattering). The situation is reversed 
for the valley-flipping cross sections, a T f, which display a 
maximum at 0 = 7r, and vanish at 0 = 0. It is then clear 
that any backscattering in the system is caused solely by 
valley flipping events. 

The scattering matrix elements are naturally related to 
experimentally measurable quantities via the total, trans¬ 
port and skew cross sections defined in appendix [B] The 
transport cross section a tr is related to the transport time 


Rr — 'H'imp'VF&tr 5 


(37) 


and the total cross section <j t to the elastic scattering time 
AT 1 = n irnp v F (j t , (38) 

where rq mp is the impurity concentration. The ratio 


£ _ Rr _ &t 

T e (J i r 


(39) 


reflects the degree of angular isotropy in the scattering 
processes j 35 l 36 E 5 l with £ = 1 in the fully isotropic (mas¬ 
sive) limit at low energies. The skew cross section , 
quantifies the asymmetry of the scattering about the elec¬ 
tron’s direction of incidence, and it is related to the ap¬ 
pearance of Hall currents, as characterized through the 
skew parameter 10 


7 = 


Rr 

3~sk 


&sk 

(7tr 


(40) 


The form of the differential cross section in Eq. [36] 
results in no skewness, since it is symmetric about the 
x-axis (0 = 0), yielding null skew cross sections (see ap¬ 
pendix [b| . 

The ratio of scattering times, 


£ = 


( Rpjp 

tr,rr' 


(41) 


can be determined analytically in the low energy regime 
kR < 1, where the only contributing channels to the 
scattering are the two lowest, j = Using the rela¬ 


tions in Eq. |B10[ one gets 


^ _ Rr ^ 

R 


(l/ ++ ,jl 2 +l/ + -,}l 2 ) 


I/. 


++, 


+ 3 I/+- il 


(42) 


Remote charge impurities lead to a local change in the 
carrier concentration in their area of influence and a lo¬ 
cal change of the Fermi momentum in that region. These 
impurities can then be described in the continuum model 
by the perturbation vl , which respects all the symmetries 
of the graphene lattice and belongs to the A\ representa¬ 
tion of the Cq u point group. 30 39 This shift does not lead 
to valley mixing and from Eq. 42 (as / r<fj i = 0) one finds 


£ = 2 for kR < l, 
been experimentally shown. 35 36 


regardless_ of the value of v, as has 
This is a direct conse¬ 


quence and signature of the anisotropy and predominant 
forward scattering in graphene. 

Some impurities may generate a staggered potential 
that differentiates between the A and B sites, as men¬ 
tioned above, with a perturbation described by in 
the continuum model. This term reduces the point sym¬ 
metries of the lattice (belongs to the B 2 representation 
of C§ v ), and although it does not mix valleys, it opens a 
gap. Moreover, a staggered potential does not commute 
with h or 7i, reducing the s et of constraints on the scat¬ 
tering matrix (see Eq. B12), to those from time reversal 
and unitarity. This situation leads to qualitative changes 
in the isotropy of the scattering process. The presence of 
s ^ 0 not only affects the ratio of scattering times £, bu t 
also results in a non-zero skew cross section (Eq. |B14 ), 
given that 


&sk,TT - &sk, 7 


and 


&sk,' 


f = 0 


(43) 


(44) 


and the corresponding appearance of a valley Hall 
effect.^ This effect is associated with the appearance of 
transverse valley currents, which can be characterized by 
a valley Hall angle 


&VH = 


JVH 

jx 


(45) 


At zero temperature the valley Hall angle is equal to the 
valley skewness, in the absence of side-jump effects (see 
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appendix |B|), 

lv = ^(iK ~1K>) • (46) 

This effect is found to be robust to a non-zero t valley- 
mixing perturbation. Note that in the case of valley po¬ 
larized incident currents, 7 = 7 y, leading to the appear¬ 
ance of Hall voltages due to the accumulation of charge 
on the sample edges. 

For impurities that accommodate on the center of 
the hexagon, as is the case of Ca and A1 atoms,® 
for a substrate with a periodicity commensurate with 
graphene,® and for ‘Kekule’ distortions with zero dimer¬ 
ization anglej^ the perturbation terms can be described 
by Reitfie 11 ~ t/3 ; this does not break point symme¬ 
tries of graphene (belongs to the Ai representation of 
Cq v ) but breaks translational invariance, which leads to 
vanishing Dirac points and t he open ing of a gap.® Other 
hexagon-centered impurities I 23 I 24 I 5 2 J lead to more general 
perturbations with 0 7 ^ 0 , as described by t/3e 11 ^ (in the 
B\ representation of Cq v ), which also open a gap in the 
spectrum. The presence of these interactions is reflected 
in particular through the deviation of £ away from 2 at 
low energies, as we will describe in detail below. 

In the presence of t and s perturbations, £ not only 
depends on these parameters but also on the value of 
v. Since { 0 ^ 3 , /3e zry ^} = 0 , and as c^,/3, and 7 s are self¬ 
inverse (involutory), and obey {/3,cq} = {a^a^} = 0 , 
the local spectrum shows a gap of 2 A, 

k = y/(E-v ) 2 - A 2 , (47) 


where 


It is then clear that as graphene is perturbed by 
adatoms or deformations, and different symmetry break¬ 
ing terms appear in the Hamiltonian, many of the dynam¬ 
ical electronic properties will change depending on details 
of the specific terms. In particular, the ratio £ will gener- 
ically deviate from its value of 2 as the massless nature of 
free graphene quasiparticles is affected by these pertur¬ 
bations. The appearance of backscattering becomes in¬ 
creasingly important for stronger perturbations and can 
be detected by different transport measurements. 


B. Spin-orbit interactions 

In the presence of spin-orbit interactions (SOIs) 
we consider an extended representation 0 = 

C0A,K,t ? 0 5 0A,K',t? ^B,K,i, ^B,K',i, 

^A,K',i) T ^ and define the corresponding matrices 

A = <8> So = ( , (49) 

B = f3 ® s 0 = ^ , ( 50 ) 

T 5 = 7 5 ® so = (^ i) y> ) , (51) 

which operate on the pseudospin and valley degrees of 

freedom, but leave the spin unchanged, as so is the iden¬ 
tity in spin space. These operators satisfy 

{Ai,Aj} = 2 iSij , ( 52 ) 

{B,Ai} = 0 , ( 53 ) 

[A,r 5 ] = 0 , (54) 


A 2 = s 2 + t 2 


(48) 


When the potential shift generated by the impurity is 
such that \E — v\ |A|, one recovers the linear (ultra- 
relativistic regime), k ~ \E — v\ + (D(A 2 /\E — v\), and the 
effect of A is minimal on the scattering amplitudes; this 
leads to / r ^ ± 1 « 0 , as well as to / rr i « f TT _ 1 , and 
consequently £ « 2 , even when s and t may be nonzero. 
A should then be comparable to the local kinetic energy, 
A ~ \E — v\, in order to have a noticeable effect. In this 
limit the valley-flipping contributions are non negligible, 
with / TTj 1 7 ^ f TT - 1 , and the ratio £ exhibits a large drop 
at low carrier concentrations, as can be seen in Fig. EK 
It should be noted that although the ratio £ has sim¬ 
ilar dependence on energy in the presence of s or £ the 
mechanisms by which the scattering isotropy is enhanced 
are different. In the presence of t/3e n 0 the isotropy 
is achieved through valley-flip processes that acquire an 
amplitude comparable to the valley preserving processes 
as \t\/\E-v\ 1, leading to £ « 1 in Eq. 42 In the case 


of the staggered perturbation, the isotropy is enhanced 
by the decay in amplitude of one of the two available scat¬ 
tering channels (e.g. f TT 1 0 while f TT _ 1 7 ^ 0 at the K 

point) as |s|/||E — v\ 1 , at each valley independently.® 


and 

{£,r 5 } = 0. (55) 

We also define 

Si = Si 0 (j 0 0 r 0 , (56) 

which act on the spin degree of freedom, with 

{Si, Sj} = 25ij , (57) 

and commute with Ai, B and T 5 . In this representation, 
the time reversal,^ helicity operator, and total angular 
momentum are defined as 


T = iBT 5 A 1 S 2 C , 

(58) 

h = T 5 A-k/k , 

(59) 

j z = —ihd$ + r s A/2 + S 3 / 2 , 

(60) 

with J z having integer eigenvalues, T 2 = 
{T, J z } = 0 . Similarly, one can define’ 

— 1, and 

Ti = e“ ir ^r 5 i 2 5 2 C , 

(61) 

with 77 — 1 , and satisfying 


{71,./,} = 0. 

(62) 
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FIG. 1. Numerical calculation of the transport to elastic time ratio £ (£50 in the presence of SOIs), for a set of 500 randomly 
sized impurities, where (5 < R < 15)A, and v = 2eV, for the following systems: a) In the presence of symmetry breaking 
effects described by t and s, notice £ deviates from the value of 2 by nearly a constant value at large energies, and reaches its 
minimum at/near zero energy. Inset: Dependence of £ on A = \/t 2 + s 2 for fixed values of v and E = 60meV. b) Effects of t 
and s on £50 in the presence of a Rashba interaction A r = 20meV. Notice the sharp drop of £ so at E —> 0 due to the Rashba 
interaction. The deviation from 2 at high energies is due to t and s, allowing the quantification of their amplitudes. Horizontal 
dashed lines correspond to the value £50 for given t and s in the inset of a); the vertical dashed line shows the value of A r; 
£so drops rapidly for E < Xr. c) Effects of t and s on £50 in the presence of intrinsic SOI, A so = v/10. Notice that £so has 
a similar behavior as in a), indicating that the simultaneous extraction of intrinsic SOI and the symmetry breaking amplitudes 
t and s from £so is not feasible. 


1. Symmetry considerations and constraints on spin 
transport 

In this representation, the set of spin dependent inter¬ 
actions that conserve total angular momentum is found 
to be 

V$°J = {MA 1 S 2 - A 2 S 1 ), A (A 2 S 2 + ASi), 

A so r 5 ^353,Asr 5 5 3 }, (63) 

all with real coefficients. These interactions commute 
with both T and 7i- For 



H = Ho+ V Tl , 

(64) 

where 

V Tl =V*°U{vI,tBe ir5 *} , 

(65) 

one finds 

[T t ,H} = [T, H] = 0, 

(66) 

and utilizing the unitarity of the S matrix one 
the condition 

obtains 


T^TvSjTR 1 r 5 = Slj 

(67) 


(see Eq. C5). 

It should be noted that the SO terms Xr and A do 
not lead to a gap opening; instead they break the spin 
degeneracy by doubling the bands for a fixed energy, such 
that 


( 68 ) 


where 

T = ^A^ + A* . (69) 

These interactions split the low energy resonant peaks 
corresponding to j = 0 and j = ±1 channels whenever 
T > 2E^ and produce a noticeable drop at low carrier 
concentrations (E < T) in the ratio of cross sections, 

Cso = ^rr>ss' t ^rr'.s' , ( 70 ) 

X^tt'ss' a tr,rr'ss' 

as shown in Fig. [TJd. 

The spin independent terms lead also to a variation 
of £50 at low energies, although with a different depen¬ 
dence on the carrier concentration, as shown in Fig. [l^i. 
Therefore, in the presence of Xr, A, t and s terms one 
can determine the value of v and the ratio A/v from 
the behavior of £50 with Fermi energy, as illustrated in 
Fig. [lp 

The intrinsic SOI (Aso) does not commute with T\ 
and does not cause spin flips, but opens a gap, and 
leads to a reduction of £50 5 such that £50 ^ 1 as 
Aso/\E — v\ —» 1. However, the gap generated by A so 
is incompatible with the gaps generated by t and s , since 
these interactions do not anticommute with the intrinsic 
SOlP In this case, the effect of the intrinsic SOI can 
not be distinguished from symmetry breaking effects in 
samples through measurements of the ratio £ 50 - I n the 
presence of A so and with no symmetry breaking terms 
(t = s = 0), the dependence of ^so(E) on the differ¬ 
ent values of A so are qualitatively similar to those in 
Fig. [it. It is the case that the simultaneous presence of 
A so, t and s will lead an overestimation of A so from 
£50 measurements, due to the similar deviation caused 
by t and s, as shown in Fig. [TJ^. 


k± = y/{E - v) 2 ± 2T(£i - v) 
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C. Skew scattering and spin Hall effect 


An alternative method to measure the strength of SOIs 
in graphene was introduced recently by Ferreira et al 10 
The method exploits the skew character of scattering in 
the presence of locally enhanced SOIs. Here, we study 
effects produced by dimerization perturbations such as 
tBe lT and staggered sublattice symmetry breaking, 
and how they influence skew scattering and the possi¬ 
ble determination of the SOI parameters from transport 
experiments. 

In the absence of valley mixing effects or sublattice 
symmetry breaking, the V^?j perturbations commute 
with both T and 7i, as mentioned above. The sym¬ 
metries of the scattering matrix result in cr s k,TT ft = 
&sk,T rX-l and <7sfc,TTt7 ^’s/cjTTpr 0 (see ap- 

pendixJCJ). This indicates that carriers with opposite spin 
are scattered in opposite directions, which would lead to 
their accumulation on the sample edges and the appear¬ 
ance of a spin Hall effect (SHE) signal . 10 

One way to quantify this effect is through the ratio of 
the skew to transport cross section for each spin projec¬ 
tion, 7 ^ and 77 and define the spin transport skewness 
as 


75 = 2 (7t — 7t) 


where 


7 8 = 


Err's' °sk,rr',s 

EttV &tr,TT's 1 


(71) 


(72) 


(see Eq. Cll in appendix), and 7 $ considers valley and 


spin preserving events, as well as valley and spin rever¬ 
sals. The spin skewness 7 s considers then all processes 
resulting in spin asymmetries, it is enhanced through res¬ 
onant scattering, and fully characterizes the spin Hall 
angle at zero temperature, 


= 75 , (73) 

T =0 

in the absence of side-jump effectsP^ In the absence of 
valley mixing or sublattice symmetry breaking, the skew 
parameters for each spin are equal in magnitude and op¬ 
posite in sign, leading to 7 s = 77 This effect was re¬ 
cently observed in graphene samples in which the SO 
was enhanced through the deposition of large copper or 
gold nanoparticlesp - 1 and in weakly hydrogenated sam¬ 
ples, and found to yield rather strong SHE signals ! 9 

In the presence of valley mixing, tBe lT in addition 
to SOI terms, both T and 71 commute with the Hamil¬ 
tonian, resulting in the following relations, <J s k,TT,ss = 

&sk,rr,ss — &sk,ff,ss — &sk,ff,ss cknd &sk,TT' ,SS - 

(Tsk,TT f ,ss = 0 , as well as a nonzero value of the skew 
parameter 7 ; its dependence on t as a function of carrier 
concentration is illustrated in Fig. [2^i. 

However, in the presence of a staggered perturbation 
sAs, the Hamiltonian does not commute with the opera¬ 
tor 7i • This reduces the set of conditions imposed on the 


Osh = 


Jsh 




FIG. 2. a) Spin skew parameter 7 s as function of energy for 
different values of t in a system with A so = 25meV, v — 2eV 
and R — 6 A. Even for large values of t the skew scattering 
persists, as expected from symmetry arguments, b) Map of 
spin skew parameter 7 s vs symmetry breaking staggered in¬ 
teraction s and energy (carrier concentration) E, for A so, v 
and R as in a). Notice that for large values of s the skew pa¬ 
rameter 7 —>• 0 , showing that skew scattering is not robust to 
this kind of perturbation, c) Spin skew parameter dependence 
on s for a set of 2000 randomly sized impurities (3 < R < 9) A, 
for values of v and A so as in (a). This shows the robustness 
of skew scattering to random impurity size, but not under s 
perturbations, d) Valley skew parameter dependence on s for 
the same set of random impurities and s values in (c). Note 
that in systems with s > A so, it is generally the case that 
|7v| > 1 75 1 (notice different vertical scales). Both 7 s and yy 
vanish away from the resonant regime, as it is the case for 
s/Aso = 80. 


skew cross section elements, so that the relation between 
&sk,TT,ss and (jsk , tt, ss dep ends on the staggered param¬ 
eter s (see Eqs. E3 and |Clip . A nonzero s leads to 
a reduction of 7 s and the corresponding transverse spin 
currents. Figure (2 Jd shows the energy dependence of 7 s 
for different values of the staggered perturbation s in the 
system of Fig. UK For increasing values of s the spin skew 
parameter 7 s decreases and eventually vanishes for large 
s values. In Fig. [2]: we show the effect of s on a sample 
of randomly sized impurities. We notice that the SHE is 
robust to impurity size disorder, but it is not to staggered 
potential perturbations; 75 decreases with increasing val¬ 
ues of 5 , eventually vanishing for all Fermi energies. It 
is also interesting that the presence of s leads a non-zero 
valley skewness, yy, which would then result in a valley 
Hall effect, accumulating carriers from different valleys 
to different edges of the sample. Figure |2]l shows the de¬ 
pendence of 7 v on energy and s strength, for a random 
impurity configuration, reflecting a relatively large value 
of 7 v for sizable s, while 7 s is much smaller. At high 
values of 5 , however, the resonance conditions are offset 
and both yy and 75 are small in the energy range shown. 































FIG. 3. Schematic representation of spin and valley skew scat¬ 
tering processes due to different perturbations. SOI perturba¬ 
tions act as a spin splitter while leaving the valleys unaffected. 
The skew character of the scattering remains unaffected in the 
presence of Vj\ • In contrast, the presence perturbations that 
belong to Vf, j and do not commute with 7i, such as the stag¬ 
gered potential s, result in valley skewness. In this case the 
perturbation acts as a valley splitter, while leaving the spin 
unaffected. In the presence of a staggered potential and SOIs, 
spin skewness will be reduced, depending on the strength of 
the staggered effects. 

We emphasize that SHE appears and it is enhanced 
close to the resonant scattering regimef 10 due to the pres¬ 
ence of centrosymmetric impurities that enhance SOIs 
(and produce also a local potential shift). This effect 
is robust to impurity size and potential disorder,® as 
well as to perturbations that lead to valley mixing in the 
form tBe lF ^. Notice also that close to resonance, both 
the spin and valley Hall effects are enhanced, but be¬ 
come negligible far away from resonance, as is the case 
for s/Aso = 80 in Fig- and d. The competition be¬ 
tween these two effects can be schematically illustrated 
as in Fig. [3] the SHE generated by skew scattering is 
not robust to a staggered potential, as this leads to a 
non-zero valley skewness. The latter reduces spin accu¬ 
mulation on the sample edges, while enhancing the valley 
contrast. 


III. DISCUSSION 

In light of our results and examination of symmetry 
breaking effects, let us consider recent experimental re¬ 
sults on decorated graphene. 

In the experiments by Monteverde et a/.p- 6 graphene 
was deposited on Si 02 , and a ratio of transport to elastic 
times £ ~ 1.8 was observed. The lowest carrier concen¬ 
tration for which £ is reported corresponds to E « 100 
meV (see their Fig. 3). As £ did not display a sharp drop 
below the value of 2 until the experimental minimum E 
considered, we can safely consider that the value of any 
Rashba SOI is lower than « 200 meV, if present.® [It is 
interesting to notice that sample A in their Fig. 3 exhibits 
a drop in the value of £ for both polarities near the min¬ 
imum carrier density, suggesting the presence of Rashba 


SOI in that sample.] We can also neglect the intrinsic SOI 
in pristine graphene^ 2 and it being unlikely enhanced by 
the Si02 substrate.® However, the effects of substrate 
and impurities are more likely to be symmetry breaking, 
with a staggered s, and/or hopping modulation £, char¬ 
acter. Assuming this to be the case, and comparing with 
the results of Fig. [lj), we estimate the value of A/v (with 
A = \/t 2 + s 2 ) to be A/v ~ 1/7, or A ^ 300meV - not 
an unreasonable estimate of the staggered perturbation 
strength. 

In the experiments by Marchenko et al .p Rashba in¬ 
teraction was enhanced through the intercalation of gold 
atoms between graphene and the Ni substrate. The value 
of the Rashba interaction was found to depend on the 
detailed position of the gold atoms with respect to the 
graphene sample. For samples exhibiting a Rashba inter¬ 
action of ^ 70meV, gold atoms were located in the center 
of the hexagon; for ~ 60meV, the gold atoms were in the 
center of non consecutive hexagons; and for ^ 7meV, the 
gold atoms were on top of carbon locations. Each of these 
arrangement can be seen to generate different symme¬ 
try breaking terms, in addition to the enhanced SOL In 
the case of hexagon-centered gold atoms one expects the 
perturbation to be dominated by t and u, as these terms 
possess the symmetries of this specific arrangement. In 
the cases of 60 and 7meV, one can expect the symmetry 
breaking to be described by {t,s,v} and respec¬ 

tively. In those cases, the detection and quantification of 
the Rashba SOI based on transport and elastic times not 
only allows one to quantify the strength of the Rashba in¬ 
teraction but also to obtain an estimate of the additional 
parameters of the Hamiltonian (see Fig. So- 

For example, in the case of a Rashba interaction A r = 
20meV, we notice a sharp drop in the ratio £50 for a 
Fermi energy around this value (indicated by a dashed 
vertical line) in Fig. [TJd. For values of energy larger 
than 20meV, the ratio displays a nearly constant value, 
slightly shifted down from the value of 2 (as indicated by 
horizontal dashed lines in Fig. [TJd) , with an asymptotic 
value dependent on the value of A/v. 

The decoration of graphene with heavy adatoms po¬ 
sitioned at the center of hexagons leads to the enhance¬ 
ment of intrinsic SOI.® The detection and quantification 
of this interaction through transport measurements be¬ 
comes masked by symmetry breaking effects, as shown 
in Fig. [1J. This would lead to overestimation of A so 
from measurements of £50 with varying Fermi energy. 
This suggests that different measurements, such as SHE 
should be used in conjunction in order to fully character¬ 
ize the role of the adatom perturbation. 

Skew scattering and the resulting SHE represents an 
alternative method of quantification of SOIs. The ph- 
ysisorption of nanoparticles on the surface of graphene 
may be expected to produce a potential shift v and 
enhanced SOIs with no additional symmetry breaking 
effects, given their large size when compared to the 
graphene lattice constant, leading to more accurate quan¬ 
tification of the SOIs. 11 Intrinsic SOI enhanced by heavy 
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adatoms,® may be characterized through this method. 
One can expect staggered effects to be weak in these 
samples, s <C A so, so that skew scattering would be 
an important contribution to the SHE. 

Our results are also in general agreement with the the¬ 
oretical analysis by Pachoud et a/.,® who examine the 
different low-energy perturbations arising from different 
adatom locations on the graphene lattice. They notice 
that the Hall effects produced from hexagon centered 
or on-top adatoms are of different character. For the 
hexagon centered position they expect a pure spin Hall 
effect, in agreement with the symmetry constrains pre¬ 
sented in our work. Moreover, for on-top positions it was 
suggested that if the adatom number of on -A and on -B 
types is different, a Hall effect may be observed. This is 
in agreement with our results for the staggered potential 
perturbation and the corresponding appearance of a val¬ 
ley Hall effect, which could eventually yield a charge Hall 
effect for valley polarized currents. 

Finally, in hydrogenated samples, where the Rashba 
interaction is enhanced due to sp 3 deformations,® 58 the 
non-local resistance was seen to be independent of carrier 
concentration. This suggested that skew scattering was 
negligible, and the appearance of SHE was attributed 
to the side-jump mechanism.The weak contribution of 
the skew scattering in this kind of sample, however, may 
also be due to the presence of strong symmetry breaking 
effects that reduce the conditions for the appearance of 
SHE, such as a strong staggered perturbation, s Xr. 
This regime would result in a negligible spin skew pa¬ 
rameter, as shown in Fig. [2^, due to the dominance of 
the valley skewness. The persistence of SHE would fur¬ 
ther suggest that side jump effects are robust to this kind 
of perturbations, and/or that the phenomenon observed 
should consider the role of valley skewness in the non 
local resistance. It would be interesting to carry out a 
detailed study of this issue in the future. 


IV. CONCLUSIONS 

We have studied how symmetry breaking effects lead 
to the appearance of backscattering in graphene, which 
is reflected among other things on the deviation of the 
transport to elastic time ratio £ from its ideal value of 2 
at low carrier concentrations. We have also considered 
that perturbations introduced by adatoms in graphene 
lead to the enhancement of Rashba-like interactions, in 
addition to various symmetry breaking effects. In many 
situations, we find it is possible to detect and separately 
quantify the SOI and the symmetry breaking amplitudes 
through measurements of the ratio of transport and en¬ 
ergy relaxation times, £, as function of carrier concen¬ 
tration. This separability of contributions becomes a 
challenging task if the SOI enhancement has an intrinsic 
character, however. In that case, the dependence of £ 
on carrier concentration for SOI and symmetry breaking 
interactions is quite similar, obstructing the extraction 


of each effect independently. Finally, we discussed that 
nonzero spin and/or valley skew parameters, 75 and qy, 
are expected for certain adatoms. This allows the detec¬ 
tion of both Rashba and intrinsic SOIs in the presence 
of valley mixing effects, even when the transport times 
ratio may not give definitive answers. For adatoms that 
enhance SOIs and break sublattice symmetry, we find 
a reduction of 75, which vanishes in systems where the 
staggered perturbation is larger that the SOI amplitude. 
As a consequence, the SOI strength would be masked by 
sublattice symmetry breaking, and the appearance of a 
valley Hall effect. 
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Appendix A: Spinless scattering matrix 

In the presence of the Vt,j,h interactions the oper¬ 
ator 7i = ie~' l ' y5( ^ / y 5 a 2 C and time reversal operator 
T = /37 5 aqC, where C is the complex conjugation op¬ 
erator, commute with the total Hamiltonian. In ad- 
dition, 7? = -1 , T 2 = 1, {71,./,} = { T,J Z } = 
0 , [71, ft] = 0, {T, h} = 0, which allows one to 
choose the states j and — j, such that Tifj±j = 
d=(—1 and T'ip±j = =f( — 1 Note that 

in cylindrical coordinates one can explicitly write eigen¬ 
states that explicitly satisfy these relations, such as 
i>j = (Jj-i/2,iJj+i/2e ie ,Jj-i/2,-iJj+i/2e ie ) T e iU ~^ )e , 
for E > 0, and where Ji(kr) are Bessel functions.^ 

In order to study the scattering problem, the incoming 
and outgoing states are eigenstates of the free Hamil¬ 
tonian, = E'tp. These eigenstates can be repre¬ 

sented with the help of the (helicity) projection opera¬ 
tors P± = (1 =b 7 5 )/2, so that a general eigenstate of the 
Hamiltonian can be written as 

-0 = P+fj + P-'if (Al) 

with [ft, P±\ = 0 and h(P±'if) = ± P±'ip. In addition, for 
a partial wave component we have 

Ti'ipj = TiP+ipj + TiP-'ifj 
= (-iy(p + -p-)ip. j 


= (-1) J 7 5 y~j ■ 

(A2) 

Similarly, for the operator T, one gets 


Tipj = (-1 r + 5 7 5 v>_j. 

(A3) 


10 


The scattering of a single partial wave j can be writ¬ 
ten in terms of incoming ip™ and outgoing 'ip° ut radial 
eigenstates of the free Hamiltonian, 


where, 




(A4) 


or equivalently, 

S++J 

-S- +J e- 2i * 


-S+-je 2i ^ 

S-i 


_ f S++S-+ 
S+-.-4 


,3 / \ — -3 

Similarly, symmetries under T result in 

7 5 T5 i r- 1 7 5 = si,, 


~3 

(A10) 


(All) 




5 ++,j 


S-+J s. 

Sj — S T t' ,j^o • 


+ — >i 

3 


(A5) 


Applying the operator 7i to the partial wave compo¬ 
nent ipj and using A2 we have 


Tiipj = Titf n + TiSjTT'Titf 


(A6) 


or 


^ = r-f + 7 5 riS,r 1 - 1 7 5 V’“ • 


(A7) 


Since the S' matrix takes a block diagonal form, unitarity 
SS t = 1, requires S^-Sj = 1. Comparison with the — j 
wave component, 




results in the following condition for Sj, 


7 5 TiS i r 1 - 1 7 5 = Si,, 


(A8) 


(A9) 


S--j -S- + j 

S-i— j S ++ j 

with the resulting relations 


S++,-j . j 

S+-,-i 


S'r T,j — S TT,—j — S' r7 -. —;y — S tt,.? ' 

O _ . — .p-2ir0 _ c_ . p -2iT</> _ _c _ 

k_J q-q- - k-J q-q- ^O - kJ q-q- O - kJ q-q- ^ 


(A12) 


(A13) 

“J ’ 


where 


Appendix B: Spinless cross sections 

The amplitudes and corresponding cross sections can 
be obtained by direct comparison of the far field forms 
of the incoming and scattered waves in Eq. |A4| For an 
incoming wave front along the x-axis, one can write 


^ = e ikr cos 0^. n . 


a ikr 


(Bl) 


where fj is a matrix containing the different scattering 
amplitudes for each angular momentum channel j, and 
<pi n is a spinor representing the different incoming valley 
weights. The scattering amplitudes in each j-block are 
given bjEH 


= eyy4 / ( 5 ++i . _ i) _is + - d 

Jj V^k V iS-+J (S—,i - 1) 

and related to different cross sections by 


(B2) 


<Ttt ' W_ 27 Tk 


&t,TT' 


JTT',3 | 5 


E frr'je^ 

= —oo 

/ 7T -I 

a TT ,{0)de=- Y I fr 
_7r j =-oo 

/ 7r ^ oo 

<J r r / (^)(f COS 0) d<9 &t,rr' ^ ^ ^ ? 

_7r j = ~oo 

/ 7T 1 00 

(7 TT /(6»)sin6»d6» = - E Im (frr',jfrT',j+l) , 


&sk,TT' 


3=-oo 


(B3) 
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where j is a half integer in all the sums. Here, a TT >(0) Consequently, 
is the differential cross section, while a t , r r' is the total, 

<Jtr, TT ' is the transport, and <J s k, T T' is the skew cross sec¬ 
tion, all with valley preserving (r = r') and valley flip¬ 
ping (r 7 ^ r') channels. We similarly define the valley 
averaged quantities, 


v(0) = Ett'Ctt'W , 

Et' &t,TT' 1 

@tr Et' ®tr,TT' > 

@sk ®sk,rr' • (B4) 


EjE-oo frr,je ijf> = Ej> 1 2/rrj COS (j 6 ) , (B 6 ) 

E“-oo frfje ^ 0 = Ej>i 2 ifrfj sin (jO) , 


The presence of Vt, j,h perturbations ensures that the set 
of relations in Eq. |A13| are satisfied, so that 


frr,j frr,—j frr,— j frr. 


J ? 


(B5) 


f _ . — _ f_ . P -^ T( P — r 

JTT,J JTTt—jZ' - J 1 


— 2 ir(f) _ _e_ —2 ircj) __ 


— f tt ,—j 


which results in Eq. [36] in the text. Notice also that 


Sj=-oo frT,jf*T,j+i ~ I f tt,^ 1 2 + 2 R e (f rT jf*T : j-\-i) j 

Sj = -oo frfjf*f,j+l = ~ I/tt,|| 2 + Sj>l ^R e {frr,j f*f ,j + l) > 


(B7) 


which results in 

e~- 


)= 0 . 


00 

J=-oo 

Ej = -oo frr,j-\-l 


(B 8 ) 


Now, when /ci? <C 1, the only scattering channels with 

1 

2 ’ 


non vanishing amplitudes are j = ±|, and the ratio of 
transport to elastic times becomes 


£ Efcr _ Ett' ^ ,T 

Ett' ^tr ,T 


where, 


®t,TT 

®t,TT 

&tr,TT 

&sk,TT 


= 2(7, 


§l'" - 

&sk,TT 0 * 


tr,TT 5 


This yields 




(i/ ++ , s p+i/ + , s p) 


°tr |/++,i | 2 + 3|/ + _i| 
which appears as Eq. 02] in the text. 


(B9) 


(BIO) 


(Bll) 


1. Vt,j perturbations 

In the presence of interactions that belong to the set 
Vt, j, the conditions on the scattering matrix elements are 


those resulting only from time reversal symmetry, such 
that they are reduced to those in |A 12 [ 

(B12) 


J TT,J 


O 

^TT ,—J •) 


Rrr,j 5 

^Tf,j Rfr,—j 1 


and correspondingly, 


frrj frr ,—j 5 

frr,j frr,—j •> 

frr,j = —frr-j • 


(B13) 


These relations cannot guarantee |B7| and |B 8 [ so that for 
example Y^=-oo Irn Urr,jfrr. J+1 ) is not identically zero 
by symmetry; instead, its value depends on the details of 
the perturbations in the Hamiltonian. As a consequence, 
one can write 


&t,TT = E,ff 1 
&tr,TT &tr,TT 5 

&sk,TT - &sk,TT 5 

E,rf E,fr 5 
&tr,rf &tr,rr 5 
&sk,TT - &sk,TT - 0 

The skew parameter is defined as 


7 = 


Err' &sk,T 


’ &tr 


(B14) 


(B15) 


The relations cr s / C)Tr 7 ^ 0 and <J s k, T T = 0 indicate the 
presence of transverse valley currents, which would result 
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in the valley Hall effect. These currents are related to the 
valley skew parameter at zero temperatures by 


where 


lv 


jvH 

jx 




r, T > &sk,r'r 

V, r > &tr,r'r 


(B16) 


(B17) 


Moreover, for valley-polarized incident flux, 7 = 7 v 4 0 
in general. This would in principle produce a charge Hall 
voltage, due to the accumulation of carriers of a given 
valley at one edge of the sample. 

This effect is absent in the usual valley unpolarized 
currents, as (J s k,rr = —&sk,TT, ^sk,rr' — 0 and 7 = 0 , 
so that there is no net charge on the edges and no Hall 
effect is produced. 

In the case of kR <C 1, we notice from Eq. |B13| that the 
amplitudes of the dominant two j-channels in this regime 
(j = ±|) are not all related by symmetry, which results 
in 1 < £ < 2 , and dependent on the specific parameters of 
the Hamiltonian. When A/\E — v\ —>■ 0 (ultrarelativistic 
regime), one recovers £ —)> 2, while £ -7- 1 as A/\E — v\ 

1 .1451551 


Appendix C: Scattering matrix and spin-orbit 
interactions 


In the presence of spin-orbit interactions that com¬ 
mute with the total angular momentum, Vj? j, and spin 
inde pendent interactions that belong to Vq\ (see Sec. 
HB), both Ti = e- iT ^Y b A 2 S 2 C and T = iBY b AiS 2 C 
operators commute with H = Hq + V 71 , and allows 
us relate the states with j and —j] here one can write 


J 


2 pj = (Jj-i , iJj , Jj- 1 e~ l ° , —iJj , Jj , iJj+ 1 e l6 , Jj , 

—iJj+le^Y'e 1 ! 9 , and where Ji(kr) are Bessel 
functions.^ 

In the presence of SO interactions, an incident state 
could be projected into valley and spin subspaces with 
the aid of the corresponding projection operators, given 
by P± = (1 ± T 5 )/2 and Pjy| ; = (1 ± S z )/ 2. A general 
eigenstate of Hq can then be written as 

4 = P+(P t 4 P^ + P-(P t + P^ . (Cl) 

The partial wave components can be chosen to sat¬ 
isfy T\P±P s i\)j = =h(-l ) J P±Ps'ip-j and TP±P s ^j = 
qF(— 1) j+ 2 P T Pg'ip_j. Consequently 

= (-1 ) j (p+(p t + p^-j - p_(p t 4 , 

= , (C2) 

and 

Ti/>j = (-1 ) j+1 HP+(Pt + PM-i - P-(Pt + Pi))1>-j , 

= (-ly+tr 6 ^. (C3) 

As before, the incoming and outgoing partial waves 

describing the scattering problem can be written as 

^^f 4 s^r . (C4) 

In this case, the scattering matrix can be written in terms 
of 4 x 4, Sj-blocks, with elements S TT ' iSS 'j. By applying 
the operator 7i to Eq. |C4[ we get 

r 5 Ti SjT^r 5 = sV, (C5) 

and by applying T we get 

T b TSjT- l T h = S ] _ j . (C6) 

From these conditions and unitarity of Sj we have, 


Brr,ss,j $TT,ss,—j ^rr,ss,—j 
Q _ . _ _C_ . _ _C_ __ ir<f> 

^TT,SS,J ^TT,SS,—J ^TTiSS,-]*-' 

Q _ _ . — C _ _ . — _Q_ _ . — _Q_ _ _ 

^rr,ss,— j ^rr,ss,— j ^tt,ss,— j c 

Q _ . — Q _ . — _Q _ 

^rr,ss,j — ^rr,ss, —j — ^tt,ss,— 


_C _ . p 2iT(£ 

^TT,SS,— J C 


(C7) 


The scattering amplitudes are given by 


m = 


g— 27r/4 

\[2i rk 


( s ++n,j ~!) ~ iS +-,n,j ~i s ++,n,j ~ s +-,uj \ 

7—KTT.i (‘S'——tt,i 1) ^ —M-U) — 7—,1\l -,j 

iS ++MJ S +-,it,j ( S ++UJ ~ !) ~ iS +-,U,j 

\ ~ S -+MJ iS —,4-tJ iS -+,U,j ( S — U,j ~ !) / 


(C 8 ) 
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and the different cross sections 




53 fj’ TT '’ s 




ij6 


\j=-oo 

oo 


(C9) 


&t,rr',ss' — ^ ^ \fj,rr',ss'\ •> 

j=~oo 


1 00 ^ 

&tr,TT f ,ss' &t,TT f ,ss' T ^ ^ f j-\~1 tt' ,SS' ) ? 


&sk 


^ oo 


I 

These help define the ratio of transport to elastic times For interactions commuting with T and 7i, we have 
in the presence of SOI, the relations 


€so = 


and skew parameter 


^rr' ,ss' 


Srr'jss' a tr,TT' ,ss' 


7 5 = 2 (7t — 7t) 


(CIO) 


(Cll) 


7* 


’y ->T ; T / ; g / &sk,' 


y, r r ' s ' &tr,TT f ,s f s 


fTT,SS,j f'TT,SS,j-\-l ~ ( fTT ,SS,j fTT ,SS ,j ~\~l) ~ ( fTT,SS,j f T T,SS,j + l) ~ /rr ,SS,J /ff ,ssj + l 5 

j — — oo j= — oo j= — OO j= — oo 

oo oo oo oo 

yy f tt, ss,jfrf,ss,j-\-l = yy frr,ssjffT,ss,j -(4 = ( y^ fTT,ss,jf tt,ss,j + i) = ( E /tt,ssj/tt,ssJ + i) 5 
j= — oo j= — oo j= — OO j= — oo 

oo oo 

E TT> i s ^ij f tt' ,SS,j-\-l) = E Tr ' :^S,j f TT f ,SS,j-\-l) = 0 5 (Cl 2 ) 


J = -oo 


J = -oo 


and consequently, 


&t,TT, 

ss 

- &t,TT,SS - 

&t,TT,SS ~ 

■ &t,TT,SS 1 (C13) 

&tr,TT, 

ss 

- &tr,TT,SS — 

■ &tr,TT,SS 

— &tr,ff,ss 5 

&sk,TT. 

ss 

&sk,TT,SS 

&sk,TT,SS &sk,TT,SS 5 

&t,TT. 

ss 

: &t,TT,SS = 

&t,TT,SS = 

: &t,TT,SS 1 

&tr,TT. 

ss 

&tr,TT,SS — 

' &tr,fr,ss 

&tr,rf ,SS 1 

&sk,TT. 

ss 

&sk,TT,SS - 

~ &sk,TT. 

,SS &sk,TT,SS 5 

&t,rr'. 

ss 

- &t,TT',SS 5 



&tr,TT f . 

ss 

&tr,TT' ,ss •> 



&sk,TT' . 

ss 

= &sk,TT' ,SS 

= 0 . 



accumulation of different spins on edges of the sample, 
and the corresponding appearance of the spin Hall effect. 


In the presence of SOIs the total angular momentum 
is an integer (j E Z), as noted before. Then for <C 1, 
there are £/^ree dominant angular momentum channels 
in the scattering process, j = 0, ±1. In the case of reso¬ 
nance, these channels resonate at different energies, lead¬ 
ing to a drop of £so towards 1 at these energies, as the 
scattering process is predominantly caused by only one 
of the channels for a given spin kind. 


This leads to 7 ^ = — 7 ^ and 7 s = 7 ^, so that even 

for a valley and spin unpolarized incident wave, there is In the case where the perturbations present do not 
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commute with 7i, instead of |C13| we have 


&t,TT,SS &t,TT,SS 1 

&tr,rr,ss — &tr,TT,SS 5 

&sk,TT,SS &sk,TT,SS 1 

&t,TT,SS = &t,TT,SS 1 
&tr,TT,SS = &tr,rf,ss •> 
&sk,TT ,SS &sk,TT,SS 1 

&sk,rr' ,ss — &sk,rr' ,ss — 0 • 


In this case, 7 s gets reduced from its optimal value, 
which indicates that Vj- j perturbations that do no com- 

(C14) 

mute with T\ lead to a reduction in the possible SHE: the 
competition with valley skew scattering (produced by a 
staggered perturbation, for example) will result in a re¬ 
duction of the spin imbalance accumulated at the edges 
of the samples. As such Vj-j perturbations increase, 7 $ 
may even become negligible, as 7 ^ ~ 7 ^ if valley skewness 
dominates, as schematically represented in Fig. [3j 
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